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What atomic liquids can teach us about quark liquids 

Thomas SCHAFER*) 
Department of Physics, North Carolina State University, Raleigh, NC 27695 

", We discuss some aspects of cold atomic Fermi gases in the unitarity limit that are of 

interest in connection with the physics of quark matter and the quark gluon plasma. We 
' consider, in particular, the equation of state, transport properties, the critical temperature 

| for pair condensation, and the response to a pair breaking stress. 

£ ■ 

' §1. Introduction 

CO \ Over the last ten years there has been truly remarkable progress in the study of 

cold, dilute gases of fermionic atoms in which the scattering length a of the atoms 
__i . can be controlled experimentally. These systems can be realized in the laboratory 

^ \ using Feshbach resonances, see 1 ) for a review. A small negative scattering length 

corresponds to a weak attractive interaction between the atoms. This case is known 
as the BCS (Bardeen-Cooper-Schrieffer) limit. As the strength of the interaction 
increases the scattering length becomes larger. It diverges at the point where a 
bound state is formed. The point a = oo is called the unitarity limit, because the 
scattering cross section saturates the s-wave unitarity bound a = Air/k 2 . On the 
other side of the resonance the scattering length is positive. In the BEC (Bose- 
Einstein condensation) limit the interaction is strongly attractive and the fermions 
^ ■ form deeply bound molecules. 

<]} [ The unitarity limit is of particular interest. In this limit the atoms form a 

i-^H strongly coupled quantum liquid which exhibits universal behavior. In the BCS 

. limit the atomic gas is characterized by the small parameter (kpa), where kp is the 

Fermi momentum. In the unitarity limit this parameter is infinite and the system 
is strongly coupled. Universality arises from the fact that short distance effects are 
suppressed by the small parameter (kpr), where r is the effective range. In the 
experiments performed to date (kpr) <C 1, and we expect that the theoretical limit 
(hp a) — ► oo, (kpr) — ► is well defined. Dilute fermions in the unitarity limit provide 
an interesting model system in which a number of questions regarding the behavior 
of strongly coupled quantum liquids can be studied. In this contribution we shall 
study a number of issues that are of relevance to the phase diagram of QCD: 

• What is the equation of state at strong coupling? Is the transition from weak 
to strong coupling smooth? 

• What are the transport properties at strong coupling? Are transport properties 
more sensitive to the coupling than thermodynamic quantities? Do atomic 
liquids respect the proposed bound on the ratio of shear viscosity to entropy 
density? 
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• What is the critical temperature for pairing? Is there a universal upper bound 
on T c /Ep, where T c is the critical temperature and Ep is the Fermi energy. 

• How does the paired state below T c respond to a pair breaking stress? Are there 
any intermediate states that separate the fully paired state from the normal 
state? 

§2. Equation of State 

Asymptotic freedom implies that the equation of state of a quark gluon plasma 
at T 3> Aqcd is that of a free gas of quarks and gluons. Numerical results from 
lattice QCD calculations show that at T ~ 2T C , which is relevant to the early stages 
of heavy ion collisions at RHIC, the pressure and energy density reach about 85% 
of the free gas limit. This is consistent with the first order perturbative correction. 
Higher order terms in the perturbative expansion are very poorly convergent, but 
this problem can be addressed using resummation techniques. 2 ) In this framework 
the degrees of freedom are dressed quasi-quarks and quasi-gluons, and these quasi- 
particles are weakly interacting. 

Transport properties of the plasma indicate that this may not be correct. Ex- 
periments at RHIC indicate that the viscosity of the plasma is very small, and that 
the opacity for high energy jets is very large. An interesting perspective on this 
issue is provided by a strong coupling calculation performed in the large N c limit of 
M = 4 SUSY Yang Mills theory. The calculation is based on the duality between the 
strongly coupled gauge theory and weakly coupled string theory on AdS§ x S$ dis- 
covered by Maldacena. 3 ) The correspondence can be extended to finite temperature. 
In this case the relevant configurations is an AdS^, black hole. The temperature of 
the gauge theory is given by the Hawking temperature of the black hole, and the 
entropy is given by the Hawking-Beckenstein formula. The result is that the entropy 
density of the strongly coupled field theory is equal to 3/4 of the free field theory 
value. 4 ) This implies that thermodynamics is not drastically effected in going from 
weak to strong coupling. 

The crossover from weak to strong coupling can also be studied in the context 
of cold atomic gases. In the non- interacting limit the energy per particle is given 
by E/N = 3Ep/5. The Fermi energy Ep = kp/(2m) is related to the density 
N/V = k F /(3n 2 ). In the BCS limit interactions reduce the energy per particle. To 
leading order in (kpa) we have 

E 3E F f 10 . ] 

In the unitarity limit (kpa) — > oo the energy per particle must be a universal con- 
stant times the free Fermi gas result, E/N = £(3Ep/5). The calculation of the 
dimensionless quantity £ is a non-perturbative problem. The most accurate results 
for £ are believed to come from Green Function Monte Carlo (GFMC) calculations. 
Carlson et al. find 5 ) £ = 0.44. This value is consistent with recent experimental de- 
terminations. GFMC calculations also show that the crossover from weak to strong 
coupling is smooth. 
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It is interesting to find analytical approaches to the equation of state in the 
unitarity limit. We have recently summarized this subject in. 6 ) There are a number 
of methods that can be systematically improved: 

• An expansion in the number of species. 7 ) At leading order this is essentially the 
BCS approximation, and higher orders take into account fluctuations around 
the BCS mean field. The result is £ = 0.59 + 0(l/N). 

• An epsilon expansion 8 )' 9 - 1 around d = 4 or d = 2. The 4 — e expansion involves 
weakly coupled bosons and fermions, while the 2 + e expansion is related to 
the perturbative (kpa) expansion. The most reliable results are obtained by 
combining the two methods. Arnold 10 ) et al. conclude that £ = (0.30 — 0.37). 

• An expansion in one over the number of spatial dimensions. 11 )' 12 ) This method 
corresponds to the hole line expansion of Bethe and Brueckner. The pair con- 
densation energy is formally suppressed by 1/d. The leading order result is 
f = l/2 + 0(l/d). 

The various methods emphasize different aspects of the physics of a cold fermion 
gas, and they all have their advantages and disadvantages. None of them appear to 
converge rapidly. 



§3. Transport Properties 



The matter produced at RHIC is characterized by strong radial and elliptic 
flow. 13 ) This observation has lead to the conclusion that the shear viscosity to 
entropy density ratio of the quark gluon plasma at temperatures near T c must be very 
small, 14 ) rj/s <C 1. In the weak coupling limit the shear viscosity can be computed 
in perturbative QCD. The result is 15 ) 

a = 

8 5 4 log(2.42 5 - 1 )' 1 > 

In a weakly coupled (g ~ 1) QGP rj/s is very large. The shear viscosity to entropy 
density ratio decreases as the coupling increases, but it is hard to extrapolate the 
weak coupling result to the strong coupling domain. Kovtun et al. conjectured that 
there is a universal lower bound 16 ) rj/s > ^/(4ttA;_b). The bound is saturated in the 
case of strongly coupled gauge theories that have a gravity dual, like the M = 4 
SUSY gauge theory discussed in Sect. [2 

In cold atomic gases we can reliably compute rj/s in the BCS limit. The ratio 
is temperature dependent and has a minimum at T ~ Tp, where Tp = Ep/ks 
is the Fermi temperature. The shear viscosity is proportional to 1/a 2 , and rj/s is 
very large in the weak coupling limit. As in the QCD case, there are no controlled 
strong coupling calculations. It is possible, however, to reliably extract r]/s from 
experimental data on the damping of collective oscillations. 17 ) Collective modes have 
been studied in a number of experiments. 18 ) In the unitarity limit the frequency of 
collective modes is well described by ideal hydrodynamics. 19 ) The energy dissipated 
due to viscous effects is 

E = ~ J d 3 xr/ (diVj + djVi - ^5ijd k vA - j d 3 x ((diVi) 2 , (3-2) 
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Fig. 1. Viscosity to entropy density ratio of a cold atomic gas in the unitarity limit. This plot is 
based on the damping data published in 20 ) and the thermodynamic data in. 21 )' 22 - 1 The dashed 
line shows the conjectured viscosity bound rj/s = l/(47r). 



where V{ is the flow velocity, rj is the shear viscosity and £ is the bulk viscosity. In 
the unitarity limit the system is scale invariant and the bulk viscosity in the normal 
phase vanishes. 

We recently analyzed the experimental data of the Duke group. 23 -* Kinast et 
al. measure the damping rate of the radial breathing mode in an elongated, axially 
symmetric trap. 20 ) They report measurements of the damping rate T in units of the 
radial trap frequency as a function of T/Tp. The shear viscosity to entropy density 
ratio is given by 

where ui = (uj^uJz) 1 ^ is the geometric mean of the trap frequencies, N is the number 
of atoms, and S/N is the entropy per particle. Our results are shown in Fig. [TJ Wc 
observe that rj/s has a shallow minimum near T c ~ Tp/3. The value at the minimum 
is rj/s ~ 1/3, roughly four times bigger than the proposed bound. 



§4. Critical Temperature 



In the limit of high baryon density and low temperature the QCD phase dia- 
gram contains a number of color superconducting phases. Color superconductivity 
is characterized by the formation of quark Cooper pairs. At asymptotically large 
density the attraction is due to one-gluon exchange. In this limit the pairing gap is 
given 

by 24 ' 

where g is the running coupling constant evaluated at the scale /j, and Abcs = 
2567r 4 (2/iV / ) 5 /W Here, /i is the baryon chemical potential and Nf is the number 
of flavors. This result exhibits a non-BCS like dependence on the coupling constant 
which is related to the presence of unscreened magnetic gluon exchanges. The critical 
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temperature is nevertheless given by the BCS result T c = e^A/ir. 

In the weak coupling limit the gap and the critical temperature are exponentially 
small. The ratio T c jEp increases with g and reaches a maximum of T c = 0.025Ep 
at g = 4.2. The maximum occurs at strong coupling and the result is not reliable. 
Using phenomenological interactions, or extrapolating the QCD Dyson- Schwinger 
equations into the strong coupling domain, 25 ) one finds critical temperatures as large 
as T c = 0.15£f. 

At low temperature the atomic gas becomes superfluid. If the coupling is weak 
then the gap and the critical temperature can be calculated using BCS theory. The 
result is 

^ = p|^ exp (-^)' <4 ' 2) 

where the factor (4e) 1//3 is the screening correction first computed by Gorkov et al. 26 ) 
Higher order corrections are suppressed by powers of (kpa). In BCS theory the 
critical temperature is given by T c = e^A/ir. Clearly, the critical temperature grows 
with the scattering length. Naively extrapolating equ. (|4-2p to the unitarity limit 
gives T c ~ 0.28Ep. The value of T c has been determined in a number of Monte Carlo 
calculations. Burovski et al. find 27 ) T c = 0.152(7)£'f, while Bulgac et al. obtain 28 ) 
T c = 0.23(2)E F and Akkinei et al. quote 29 ) T c = 0.25E F . The larger values of T c are 
in better agreement with the transition observed in trapped systems. 30 ) 

§5. Stressed Pairing 

The exact nature of the color superconducting phase in QCD depends on the 
baryon chemical potential, the number of quark flavors and on their masses. If the 
baryon chemical is much larger than the quark masses then the ground state of 
QCD with three flavors is the color-flavor-locked (CFL) phase. The CFL phase is 
characterized by the pair condensate 31 ) 

d>?C l5 rf) = (6?6 b j -6pl)<t>. (5-1) 

This condensate leads to a gap in the excitation spectrum of all fermions and com- 
pletely screens the gluonic interaction. Both the chiral SU(3)l x SU(3)r and color 
SU(3) symmetry are broken, but a vector-like SU(3) flavor symmetry remains un- 
broken. In the real world the strange quark mass is not equal to the masses of the up 
and down quark and flavor symmetry is broken. At high baryon density the effect of 
the strange quark mass is governed by the shift fi s = m 2 /(2/i) of the strange quark 
Fermi energy. 

The response of the CFL state to a non-zero /i s is a difficult problem that has 
not been fully resolved, even in the weak coupling limit. There are three energy 
scales that are important 

• rrix ~ (m u m s ) 1 ^ 2 (A/fi) <C A is the mass of the neutral strange Goldstone 
boson, the K . When /i s > mx the CFL phase undergoes a transition to a 
phase with kaon condensation. 32 ) 

A is the critical value of fi s at which the first fermion mode becomes 
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Fig. 2. Conjectured phase diagram for a polarized cold atomic Fermi g £LS EIS 8, function of the 
scattering length a and the difference in the chemical potentials 8{i = fj,f — /ij., from Son & 
Stephanov (2005). 



gap less. For fi s > fig the CFL phase (with or without kaon condensation) is a 



• fig ~ 2Z\ is the critical value of fi s beyond which CFL pairing breaks down 

(2) 

completely. For fi s > fi s the CFL phase is replaced by a less symmetric phase, 
like the 2SC phase or single flavor pairing in the spin-one channel. 
The most difficult part of the phase diagram is the region fii 1 ^ < fi s < fip . 
Gapless fermion modes cause instabilities in the superfluid density and the magnetic 
screening masses. 34 -* Near fig this instability can be resolved by a small Goldstone 
boson current. 35 ) Closer to fig the Goldstone boson current may become large, and 
multiple currents can appear. In this limit the ground state is more appropriately 
described as a LOFF phase. 36 ) We shall describe the LOFF state in more detail 
below. 

The atomic superfluid involves equal numbers of spin up and spin down fermions. 
The physical situation which is analogous to the response of the CFL phase to m s 
is the response of the atomic superfluid to a non-zero chemical potential coupled 
to the third component of spin, 8fi = /x-f — A conjectured (and, most likely, 
oversimplified) phase diagram for a polarized gas is shown in Fig. [2j In the BEC 
limit the gas consists of tightly bound spin singlet molecules. Adding an extra up 
or down spin requires energy A. For \5fi\ > A the system is a homogeneous mixture 
of a Bose condensate and a fully polarized Fermi gas. This mixture is stable with 
respect to phase separation. 

In the BCS limit the problem was first analyzed by Larkin, Ovchninikov, Fulde 
and Ferell (LOFF). 37 ) Consider the homogeneous solutions to the BCS gap equation 
for 5fi ^ 0. In the regime 5fi < Aq where Aq = A(Sfi = 0) the gap equation has 
a solution with gap parameter A = Aq. This solution is stable if dfi < Ao/y/2 but 
only meta-stable in the regime Aq/\/2 < Sfi < Aq. The BCS solution has vanishing 
polarization. The transition to a polarized normal phase is first order, and systems 
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at intermediate polarization correspond to mixed phases. 

LOFF studied whether it is possible to find a stable solution in which the gap 
has a spatially varying phase 

A{x) = Ae 2i ^. (5-2) 

This solution exists in the LOFF window 5/j.i < 5/j, < 5^2 with 5^i\ = Aq/\/2 ~ 
0.71Z\o and <5/i2 — 0.754Z\o- The LOFF momentum q depends on 5fi. Near 5fj,2 we 
have qvF — 1.2c)/x, where vp is the Fermi velocity. The gap A goes to zero near b\i2 
and reaches A ~ 0.25Z\o at Sfj,\. 

These results suggest that for some value of the scattering length between the 
BEC and BCS limits the homogeneous superfluid becomes unstable with respect to 
the formation of a non-zero supercurrent Vy, where (p is the phase of the condensate. 
We can study the onset of the instability using the effective lagrangian 

C = (id - e(-id) - i0<p) • A) ^ + - £(^) 2 . (5.3) 

Here, if) describes a fermion with dispersion law e(p) and p> is the superfluid Goldstone 
mode. The low energy parameters ft and / are related to the density and the velocity 
of sound. Setting up a current v s = d<p/m requires energy f 2 m 2 v 2 /2. The fermion 
dispersion law in the presence of a non-zero current is e v (p) = e(p) + v s -p — 6fj,, and a 
current can lower the energy of the fermions. The free energy functional was analyzed 
by Son and Stephanov. 38 ) They noticed that the stability of the homogeneous phase 
depends crucially on the nature of the dispersion law e(p). In the BEC limit the 
minimum of the dispersion curve is at p = and there is no current instability. In 
the BCS limit the minimum is at p ^ 0, the fermion contribution is amplified by a 
finite density of states on the Fermi surface, and the system is unstable with respect 
to the formation of a non-zero current. 

There is an ongoing effort dedicated to the experimental study of the phase 
diagram as a function of scattering length, polarization, and temperature. We cannot 
adequately summarize all of these experiments here. The MIT group has mapped out 
the superfluid-normal transition line in the (<5yu, a) plane. 39 ) In the unitarity limit the 
transition occurs at a population imbalance (nt — n^)/(n^ + n\) ~ 70 %. Currently, 
there is no evidence for inhomogeneous or anisotropic states like the LOFF phase. 

Acknowledgments: This work is supported in part by the US Department of 
Energy grant DE-FG02-03ER41260. 
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